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Abstract. In this paper we shall define the renormalization of the multiple g-zeta values (MgZV) 
which are special values of multiple q-zeta functions (q( s i> • • • » s d) when the arguments are all 
positive integers or all non-positive integers. This generalizes the work of Guo and Zhang [12) on 
the renormalization of Euler-Zagier multiple zeta values. We show that our renormalization process 
produces the same values if the MgZVs are well-defined originally and that these renormalizations of 
MqZV satisfy the q-stuffle relations if we use shifted-renormalizations for all divergent Cq( s li ■ ■• i s d) 
(i.e., S\ < 1). Moreover, when q f 1 our renormalizations agree with those of Guo and Zhang. 
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1 Introduction 



The Euler-Zagier multiple zeta functions are defined as nested generalizations of the Riemann zeta 
function: 

C( Sl ,...,s d ):= k i Sl --- k d Sd (1) 

fei>->fcd>0 



for complex variables si 



, Sd satisfying 3?(si 



+ s 



) > j for all j = 1, . . . , d. The special 



values of this function at positive integers are called multiple zeta values (MZVs) and were first 
studied systematically by Euler [9] when d = 2. Nevertheless, only in the past fifteen years these 
values have been found to have significant arithmetic, algebraic, geometric and physics meanings 
and have since been under intensive investigation (see [51 [TTJ1 fT5l fT5l |2"T] ) . 

In another direction in [22] we show by using generalized functions that multiple zeta functions 
can be analytically continued to C d as a meromorphic function with simple poles. We will henceforth 
always refer to this analytic continuation when we speak of multiple zeta functions in the rest of 
this paper. The precise location of the simple poles form the following set (see [3]): 



(si,. 



,Sd) 6 



Si = 1; , or si + s 2 £ {1} U 2Z<i, or 



€ Z<j for 3 < j < d 



(2) 



Hence MZVs at non-positive integers are not always defined. Recently, Guo and his collaborators 
([El [H]) have applied the Rota-Baxter algebra technique to the study of MZVs after noticing 
that the stuffle (stuffing+shufHe) relations reflect exactly the Rota-Baxter property. In [12] the 
renormalization is carried out for the MZVs and they show that when £(si, . . . , Sd) is defined then 
its renormalization agrees with the value itself, provided that Sj's are all positive or all non-positive. 
Moreover, these renormalizations satisfy the stuffle (or quasi-shuffle) relations. The importance of 
this result is related to the conjecture (see [15]) that to obtain all the relations among MZVs of the 
same weight it suffices to use all the double shuffle relations including those of the renormalization 
of MZVs at positive integers. 

On the other hand, we can define the g-analog (0 < q < 1) of multiple zeta functions as follows 
(see [23])- For complex variables Si, . . . , Sd satisfying 3?(si H h Sj) > j for all j = 1, . . . , d, set 



Cq(su ■■.,Sd):= 



E 



7 fel(si — 1)H \-k d (s d -l) 



ki>—>k d >0 



[k d ] Sd 



(3) 
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where for any real number r we write [r] = (1 — — q). When d = 1 this is the same as the 

g-analog of the Riemann zeta function defined in [16] . By using Euler-Maclaurin summations we 
obtained their meromorphic continuations to C d with following singularities (which are all simple) : 



(s u ...,s d ) e C d 



27T2 27T2 
■Si E 1 + : — Z, or Si E Z< + : — Zj , 

lng ~~ Ing 

27TZ 

or Si + ■ • • + s,- G Z<„- H Z for j > 1 

lng 



(4) 



Here the last part in & d is vacuous if d = 1. One can see that these g-analogues have much more 
poles than their ordinary counterparts. But when q approaches 1 we indeed recover exactly the 
poles of the multiple zeta functions. In fact, by [23] Main Theorem] for all [s\, . . . , Sd) £ C d \ &d 
lim g |i Cj( s ij . . • , Sd) — C( s ij ■ ■ ■ i s d), which shows that our g-analogue is the correct choice. 

The analytical continuation of multiple zeta functions [22j utilizes generalized functions. But 
Euler-Maclaurin summation can also be used instead which actually provides the main idea of 
special value computations contained in this paper. For future reference we define the Bernoulli 
polynomials B k {x) and its periodic analogue Bk(x) by 



te 1 



OO 

E 

n=0 



B n (x)— } , B k {x)=B k {{x}),x>l, 



(5) 



where {x} is the fractional part of x. We can then prove the analytic continuation of C( s ii • • • > s d) 
using these functions. See [23l Theorem 3.2] for more details. 

Similarly, the analytic continuation of multiple q-zeta functions Cg( s i> . . . , s d ) can be obtained 
by using Euler-Maclaurin summation formula. The major difference between ordinary MZVs and 
MqZVs is the appearance of the shifting operators Sj (1 < j < d) defined as follows: 



§>j(q(si, ...,S d ) := Cq(si, ...,S d ) + (1 - q)(q(si, ■ ■ ■ , Sj 

In general we may iterate the operator and get 



l,...,s d ). 



r, . . .,s d ). 



(6) 



Using these operators we proved the analytic continuation of multiple q-zeta functions in 23] . 

The g-analogue of MZVs will be called multiple q-zeta values (MqZVs). In this paper we will 
consider the renormalization problem for MqZVs motivated by the ideas of Guo and Zhang in 
|12j . From physics point of view these values can be regarded as the quantumization of MZVs. 
Furthermore these MqZVs also have number theoretical significance. For instance, it is well known 
that ((0) = —1/2 and £(1 — 2n) — — £?2n/(2n) for positive integers n where i?2n are Bernoulli 
numbers defined by xj (e x — 1) = J2iLo BiX l /i\. What are the right g-analogue of these numbers? It 
turns out that one of the ways to find the answer is to consider Rieman qZVs at negative integers (see 
[TBI (6)]), which shows that the odd indexed g-analogues of Bernoulli numbers are actually nonzero. 
Is it possible to generalize Bernoulli numbers to multi-Bernoulli numbers and their g-analogues? 
Maybe this problem can be solved when we carry out further studies of the renormalization of 
MZVs at negative integers. 

The major behavioral difference between MZVs and MgZVs is the appearance of the shifting 
operator in the g-analogues defined by ^ . As we mentioned in the above it is very fruitful to study 
the stuffle relations between MZVs. The g-analogue of this is a little more complicated because of 
the shifting operator which can still be handled by setting things up carefully. The main result of 
this paper is that we can define the renormalization of MqZVs when s^'s are all positive integers 
or all non-positive integers such that (i) they coincide with the MqZV if it is defined originally, 
(ii) they satisfy a shifted version of g-stufHe relations, and (iii) they become the renormalization of 
MZVs defined by Guo and Zhang in [12] when q f 1. 

To conclude this introduction we remark that currently there are two ways to order the variables 
in MVZs and our £( s ii • • • j s d) in this paper is denoted by ((s d , ■ ■ ■ , s\) in |23j . We change our 
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notation system because it is more convenient now for the readers to compare results in this paper 
to their classical counterparts in [12] which serves as the major motivation for us. 

I would like to thank Li Guo and Bin Zhang for their interest in this work and many valuable 
comments of the first draft of this paper. Thanks are also due to the referees for their many detailed 
suggestions which have greatly improved the exposition of this paper. 



2 The Rota-Baxter algebra and the g-stuffle product 



Let k be a subring of C which is usually taken to be K or C. For any fixed A 6 k a Rota-Baxter 
k-algebra of weight A (previously called a Baxter algebra) is a pair (R, P) in which R is a k-algebra 
and P : R — > R is a k-linear map, such that 



P{x)P{y) = P(xP(y)) + P{P(x)y) + XP{xy), Vx, y e R. 



(7) 



In this paper we are going to concentrate on the following two examples, both of which are 
contained in [12] . 

Example 2.1. Let £ be a complex variable such that 5R(e) < 0. Let C{{e,e -1 } be the algebra of 
convergent Laurent series in a neighborhood of e = with at worst finite order pole at 0. Write 
T = — ln(— e) which is transcendental over C{{e,e -1 }. Then we can regard R :— C{{e,e _1 }[T] 
as the polynomial algebra with the variable T and with coefficients in C{{e,e -1 }. Let P be the 
operator on R which takes the pole part. Then it's not hard to verify that (R, P) is a Rota-Baxter 
C- algebra of weight — 1. 

Example 2.2. Let !K be a connected filtered Hopf algebra over k (see [12] §2.1] for the definition) 
and let (R, P) be a commutative Rota-Baxter algebra of weight A. Define the k-algebra 3? := 
Homk(JC, R) of linear maps from J{ to R with the product compatible with the coproduct of the 
Hopf algebra JC. Then the operator CP on Hom(!K, R) defined by CP(XL) = P o L is a Rota-Baxter 
operator on 31 of weight A. This example will be used in section 2 to define the regularized MgZVs 
(see lp])L 



In the rest of this section we will construct one such Hopf algebra of Example 12.21 For any 
subset Z of C closed under addition and shifting by —1 we define the commutative semigroup 



<tt(Z):=j S n € Z> ,(s,r) E Z x M >0 | 



(8) 



with the binary operation given by 



We will only have two different choices 



for Z in this paper: Z or Z<o- The reason to require Z to be closed under shifting by —1 is because 
of the effects of shifting operators on MgZVs. To study other renormalization at other poles in the 
future we need to set Z = Z + (2ni/ \nq)Z (see (J3J). 

Define the C-bilincar pairing ( , ) on the C-algebra CVt(Z) by 



~S~ 




rs'- 




s + s - 


+ (1-9) 


■s + s' - 1" 


- T- 




Lr - 


>- 


lr + r - 


_ r + r J 



(9) 



Recall from [T2] §3.1] that we can define the algebra: 

^C<JI(Z) n 



•K 



z •= 



where 9T(Z)° = {1} is the multiplicative identity and C0 r l(Z)" is the free C-module with basis 
s yi(Z) n . Then we may equip the g-stuffle product, which is the g-analog of the quasi-shuffle product 
* for MZVs (see [H Thm. 2.1] or [H Thm. 2.2]), on 3i z as follows: for o = (a u . . . , a m ) S m(Z) m 
and b = (&i, . . . , b n ) 6 9t(Z)" we set a' = 1 if m = 1 and a' = [a^, . . . , a m ) otherwise. Then we 
a and recursively 



define l* 9 a 



a* q l 



a* q b = (ai, a'* q b) + (h, a* q b') + ((at, a'* g b') 



(10) 
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where (ai,&i) is given by ([9]). It has a connected filtered Hopf algebra structure over C when we 
define the deconcatenation coproduct suitably. If s^qS 1 — 5Z mn > (l — q) n s m ,n denotes the top 



component of 







r -/ , 


"S" 




S 








.r. 







, then wc proved in 23, Theorem 5.1] that 



(11) 



3 Regularized multiple g-zeta values 

Let's recall the classical process of renormalization. For example, let's consider the divergent series 
J^Li n, which is the series we would get if we tried to plug s = — 1 into f(s) by using definition 
([!]). We may tamper this series by multiplying a controlling factor on each term: ^n=i rie ™ e ? f° r 
some e such that SR(e) < so that we get a convergent series. By easy manipulation (see [HI (34)]) 



00 9 

E-" E = f 

n=l 



j=0 



We then call this the regularized zeta value at — 1. To recover the finite value £(— 1) we only need 
to drop the pole part 2/e 2 and then take e = 0. This process is called the renormalization. Because 
there are more than one variable in multiple zeta functions it turns out that we need to introduce a 
concept called "directional vector" (sec Definition 14. 4ft in the regularization process in order to get 
well-behaved regularized values so that the normalization works as desired. 

Turning to our MgZVs, as in the previous section we let Z be a subset of C which is closed under 
addition and shifting by —1 (which will be either Z or Z<o). For s G Z, r > 0, 3?(e) < 0, and iel 
we first define 

„ , q x ( s ~ 1 ) exp(re\x]/q x ) 
f(s,r;e,x) := pr- . 

Note that the controlling factor becomes e ex when r — 1 and q f 1. For every vector s = 



It is clear that Z„ 



([3; £ >*) 





■ ,r d ) 


"S" 








. r - 




is 


also g 



E 

n 1 >--->n d >0 j = l 



(12) 



r = (?"!,••• , r^) G (R>o) d in tH}. Following the setup of Example 12.21 we may define the C-linear 
map 



>Ri ■■= ( Z, 



S 3 X 



e <n 



( 



;e,:r) := Q^f(s d , r d ; e, f 



. :i)" z «(i 

where Q is the summation operator (denoted by P in [51 12"4"]) 

Q(/)(z) = ^/(z + n). 



"si, • 


■ , Sd-l' 








;e,x)^J 


-n, ■ 







(13) 
(14) 





S" 




( 




;s,x) 




f- 





Definition 3.1. For s G Z d and r G (M>o) d by setting x = in Z ? 

fci>— >fc d >0j=l L n 

These values are called the regularized multiple q-zeta values (at Z). 



we define 





S" 




( 




■•') 




.r- 





(15) 
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£ ] is well-defined for all s and r. In 



Because of the assumption < we see that Z q \ 

particular we don't need to restrict to a MZV-algebra as constructed in [8] §3.2] by Ebrahimi-Fard 
and Guo. Moreover, when q | 1 we recover the definition of regularized general MZV defined in 



3.1 Regularized g-Riemann zeta values 

In this subsection we deal with the g analogue of the Riemann zeta function first. Taking d = 1 
in Definition 13.11 we find that Z q (^ S ;sj = Z q (^ ^ ;re) so that it suffices to consider Z q {s;e) := 

Z q (^ 8 ; ej . We will first put Z = Z. To study these regularized values we set, similar to [IB] . 

g^^expjeg-^]) g 3 ^- 1 ) exp(£( g -* - !)/(! - g)) 

f< (x, s; e) = 7- r-r-; = ; ; 1 5t e) < 0. 

v ' (1 - g) s [x] s (1 - q x ) s K ' 

Then taking derivatives of F with respect to x we get 



(i-«r) 



aAs+2 



■exp(eg ^a:]) 



l-g V y; ^ (l-g*)^ 1 ^ V y 1 J; Vl-g/ v H> H {l-q x Y 

Let i = e{q- x - !)/(! - g), then eft = -e(lng)g- :E /(l - q)dx. We get 



(16) 



1 / e N*- 1 
^' S;£) ^ = h^ll^ 



e(8 -1)7(1-3) e i 



— di. 



(17) 



When 5 = and = 1 we have 



/oo 
F(x,0;e)efe 



q-lf-1 



lng \ £ 



q — 1 g ' 1 e ; 
lng ' -Z - 1 ' 7T 



(18) 



Simple computations yields 



F(l, 0; e) ^g- 1 exp( £ /g) = £ g^ 1 - 



i=0 

oo 



F'(l, 0; e) =(- Ing)^ 1 + £ ^V^H - 1 + 3)^r, 

f—' 1 — g t! 



;=i 



F"(x,0;e) =(lng)V* £ (V^M)' - ^(g-^])'" 1 + ^ ^ (g-[x])'- 2 )^ 



i=2 



(9-1)= 



. (] , \2„-2^3 +1 g 1 



+ (lng) 2 g x + (lng) 2 g~ 



Note that for all < g < 1 and 5R(e) < 0, Z q (0; e) = Y^nLi ^( n i 0; e) converges. By Euler-Maclaurin 
summation formula and the analytic continuation of ( q {s) given by [X 61 (12)] 



Z,(0;e) 



1 1 1 f°° - 

F(x, 0; e) dx + -F(l, 0; e) - — F'(l, 0; e) - - J B 2 (x)F"(x, 0; e) dx 



q-lf-l 
lng 



(v) + Ech) 



Z=0 



(19) 
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where B2(x) is the periodic Bernoulli polynomial defined by ([5]). Further, we have ([20l Ch. IX, 
Misc. Ex. 12]) for k > 2 



B k (x) = -k\ 



nez\{o} 



(27rm) fe 



B' k+1 (x) = (k + l)B k (x). 



(20) 



To determine the regularized normalization for Z q (s;e) for positive s we begin with the case 
8 = 1. First we have 

-1 f°° e~* 
F(x,l;e) dx = -. — / — dt. 



\nqj_c t 



Therefore 



(1-q) / F(x,l;e)dx 



q-1 f^e- 1 



In q 



dt. 

t 



(21) 



Now for any real number a > integration by parts yields 



dt = 



hit 



* In t dt = —7 — In a — 



00 1 „\l 



1=1 



(-a)' 
III 



where 7 « 0.577216 is the Euler's 7 constant. Here we have used the fact that (see pQ or [2]) 

/>OC 

T'(l) = / e-Hntdt 
Jo 



Hence 



(1-q) / F(x,l;e)dx 



i-q 

lng 



= -7- 



00 / 



Further, 



F'(l,l;e)=(lng) 



^(i.lje) =(lng) 



2 exp(e/g) 
"(l-gF)S 



(1-9) 

2 2- 3(l-g*) + (l -g*) 2 
(l-gF)3 

2 -l + 3g 



lngg — , 

1-q 1-q 

exp(eq~ x [x}) 



Consequently 

00 

Z q {l-e)=(l-q)Y j nnA\e) 
n=l 

/ r 00 1 1 1 /- 00 . 

=(l-«)(jf F(a;,l;e)da: + -i?(l,l;e)-— F'(l,l;e) - - J B 2 (x)F (x, l;e) dx 

ln(-e) + M(g) + ^7 + 0(e), 
in g In (7 



(22) 



where we have set 



A%) =9-7: 



1 g In g (1 — g)(lng)^ 



2 /-oo 



2 12g-l 



S 2 (a:) 



x\2 



2 - 3(1 - q x ) + (1 - g x ) 

(1-^)3 



■ (ice. (23) 



As a comparison we now take a look at the behavior of C, q (s) near s = 1. It is clear that 



F(x,s,0) 



qx(s-i) 
(1 - q x ) s 



G 



Thus from formula [23l (12)] we have, near 5 = 1 

oo 



n=l 



=(1-9) 
(9-1) 



1 1 1 f°° 

F(x, s, 0) dx + -F(l, s, 0) - — F'(l, s,0)-- J 



B 2 (x)F"(x,s,0)dx 



M(q) 



1 



(s — 1) hiq y ' rL/ liiq 
Taking q f 1 we have near s = 1 

CM 



■7 + 0(s - 1). 



+ limM(g) - 7 + 0(s - 1). 



1 ?Ti 



(24) 



(25) 



On the other hand, by applications of Euler-Maclaurin summation formula f |17[ p. 531] or [201 
7.21]), for all integers k, I > 1, we get 



7=1 ■' r=l K 



J3 2 z+i(x) 



,,2;+2 



Since £>2 = 1/6, putting k = I = 1, using equation (|20p and integration by parts once in (|23[) we 
get lim g |i M(q) =7 which is consistent with (f2"5)l . 
Now we can prove the following result: 



Theorem 3.2. The 



value of Z q (^ ;ej at s = 1 — n (n G Z>o ,) is 

-. \ n 00 

— nl + £&(l-n-Z) 



Z g (l - n; re) 



g-i 



T7ie vaZue 0/ .2L 



In q \rs 

; e ) at s = n E Z>o «s 



(re) 



(26) 



z=o 



Z,(n; re) = i-J^e) + (m(<z) + I— ? 7 ) . (rCi " ' 
In q \ In a / , . 

w/iere u n+ i(e) = fr(#n - M -£ )) f or n>0,H = and H n = l + \-\ \- ± for n > 1. 



(n-1) 



(re)* 



Z! 



(27) 



Proof. Equation ([26)) follows from (fl9|) by taking derivatives because -^Z q (s;e) = Z q (s — l;e). 
When s = 1 equation (|2"T|) follows from (fT9")l and (|2"2")l by integration and Abel's Theorem because 
-^Z q (l;e) = Z q (0;e). The rest follows immediately. □ 

3.2 The range of regularized MgZV 

We now turn to the general MgZVs. Although the proof of Theorem 13.31 below is similar to the 
proof of [HI Theorem 3.3] some new phenomena arise because of the shifting principle for MqZV. 
Let 



Cfte^- 1 }^ \ £ a n s n 



,n=iV 



a„eC,JVeZ 



be the algebra of Laurent series, regarded as a subalgebra of the algebra of (the germs of) mero- 
morphic functions in a neighborhood of e = with at worst finite order poles at 0. Choose lne to 
be analytic on C\(— 00, 0]. Observe that the analytic function ln(— e) on C\[0, 00) is transcendental 
over C{{£, £ -1 } by [121 Lemma 3.1] and hence there is a natural algebra injection (see [321 (19)]) 



C{{e, £-!}[- ln(-e)]— ►CpHe.e- 1 } 



(28) 



sending — ln(— e) to T. This provides an identification of C{{e,e ln(— e)] as a subalgebra of 
C[T]{{e,e-'}. 
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Theorem 3.3. (a) For any s £ (Z<o) d and r £ 





' S~ 




( 




■■') 




. f- 





(b) For any s£ Z d and f £ (R >0 ) d , Z q ( 





S~ 




( 




") 




. r - 





e) £C{{e,e- 1 }[-H-e)}. 



Proof, (a) The key is the following computation of the tail of the Z q (s; e) for s £ Z<q. First define 



exp(e[k]/q k 



k>i 



E 

fe>0 



exp(e[/c + i]/q k+l ) 



-,k+i 



Since [k + i] = [k] + q k [i] we get 



= exp(e[z]/ g ') _ exp(g[fc]/g fc+i ) = expje^V) ^ (o- £ / g <) 
y fe>0 y i 



Therefore for s £ Z< and t = — s, 

7 fc(-t-i) exp(£[A;]/g fc ) 



9 



t!/ -l^i+t-j exp(e[i]/g i ) x — « / / 
' In 9 j! 



[i]Jexp( £ M/g')C q (-Q£'" t+J " 



(29) 



Now we first prove (a) by induction on the length d. The case d = 1 corresponds to the Riemann 
g-zeta function which has been dealt with in the last section. Suppose (a) is true for length d — 1 
(e? > 1) and let s = (si, . . . , s<j). Let denote s with its i-th component removed. Let denote the 
i-th unit vector of length d — 1 with 1 at the i-th component. By definition (jT5j and the shifting 
operator ((6| for every i such that 1 < i < d, 



E 



e n 

fci>--->fed>0 m=l 

e n 



7 fc m (s m -l) 



exp(er m [fc m ]/(7 fe " 



[k>i 



,tm(«m-l) 



exp(er m [A: m ]/g fcm ) 



fei>-—>fei_l A;i_|_i>--->A: t i>0 m— X,m^i 



\h Ti 



fci_i(si-l) 



exp(er i [A: i _i]/9 fel - 1 



[k 



i-i 



+ Sjei_i 



9-1 
In g 



E 

J=0 









;z)-z q ( 


" S| - jei_i " 


'■)] 






.?i + r 4 e 4 _ 




-5 + rjei-i- 





+E 

3=0 



OC 

E 



(=-s<-j 



si - je<_i 
r; + r*e<_i 



C 9 (-0(^e)' +8<+j 



(30) 



Here if i = 1 then the terms with ej_i are 0. If i — d then the terms with ej are and the first term 
becomes Z q (^ d ;e^jZ q (^ _? By induction assumption we see that Z q y _ ; e^j £C{{e,e~ 1 }. 



Now we prove (b) by induction on length d again. 
Case 1. Si < for some 1 < i < d. Then the proof of (a) above carries over word for word here. 
Case 2. Suppose Sj > for all 1 < i < d. We use induction on the sum s := $D i=1 Si- Clearly 
s > d. If s = d, then Sj = 1 for 1 < i < d. Thus 



= E riZ « 



s - e t 
r 



;e) £C{{e,e- 1 } 



by (a). Integrating we get Z q (^ S _ £ C{{e, e ln(— e)]. The general case follows from the 
fact that C{{e, ln(— e)] is closed under integration. □ 
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Later on, for several times we are going to need the special case of ([30]) when i = 1 and i = d. 
For convenience we list them as 

Corollary 3.4. Let d > 2 be a positive integer. For any s£ (Z<o) d and f G (M>o) , 
-si e . \ | 



; g-l^(- 8l ) 
In q 



S2 -J,S 3 - 

Ti + r 2 ,r 3 , . 



,r d 



+E 

3=0 



-•si 



OO 

E § i + %( 



«2 - J, S 3 ■ 

»"i + r 2 ,r 3 , . 



, s d 



Q q {-l){r l£ ) l +^ 



(31) 



Z a 



' s~ 


;e) =Z,( 




Si, . 












. r. 






-n, . 





In g 



E 

j=0 



E 

J'=0 



J 



(s d y. / -ly-^-j 



\r d eJ 



si, . 

n, ■ 



■ , Srf-2, Sd_l - J 

,r d -2,r d -i +r d 



OO 

E 8£\*«( 



si, . 
n, ■ 



■ , Sd-2,Sd-l - J 

,r d _ 2 ,r d _i + r d 



C q (-l)(r d e) l+s ^ 

(l + Sd + j)\ 



(32) 



l=-s d —j 

The next result tells us some very useful information about the general shape of the coefficients 
of the Laurent series Z q ^ _ ; ej . This will be used crucially in the proof of the existence of the 
renormalizations of MgZVs at non-positive integers. 

Corollary 3.5. Let s G (Z<o) d and r € (M>o) d . Then the coefficient Ci of the Lauren series 



of Z q 



^""^ ge l is an M.-linear combination of rational function of the form P(r)/Q(f) 

in n,...,rd, where P, Q have no common factors. Moreover, both P and Q are homogeneous 
polynomials in n, . . . , r d , Q is the product of linear factor r\ + • • • + r d (and n if d = 2) with 
repetition allowed, and deg(P/Q) := deg(P) — deg(Q) — i. 



Proof. When d = 1 this follow from Theorem 13.21 Suppose the claim in the corollary is true when 
the length of the vector s is < d for some d > 2. It follows from (131 [) and induction assumption that 
the only possible factors in the denominator of P/Q are of the form {r\ + ■ ■ ■ + r d ) 3 and r\ for some 
j, k G Z> . But by ([32]) we see that the only possibilities are (r*i + • • • + r d ) J ', H h r<i_i) fe , and 



Hence if d > 2 then none of r\ , ri 



fd_i or r d can appear as a factor in the denominator 



of Ci. But if d = 2 then both n + T% and ri could appear (and indeed they do by the following 
formula ([35]).) 

The statements about the degrees are clear from (|30|) by induction assumption. □ 



4 Renormalization of MgZV 

Theorem 13.31 together with the map u of (|28p shows that there is an algebra homomorphism 

Z q : IK Z — CpTJfle.e- 1 } 




which restricts to an algebra homomorphism 

Z q : 5^z<o - * CftejE" 1 }. 

Note that only when there is a positive s, can .Z^ really differ from Z q . It is well-known that 
the Birkhoff decomposition (see [H] Theorem II. 5.1]) yields two maps Z q _ and Z q+ such that 
Z q = Z q _ * Z 9+ . The properties of this decomposition implies the following result immediately. 
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Proposition 4.1. The map Z q+ : 

an algebra homomorphism Z q+ : HCz <0 



-> C[T][e] is an algebra homomorphism which restricts to 
C[e]. 



To write down Z q+ explicitly we need the following definition which is slightly different from 
[H Definition 3.7]. 



Definition 4.2. Let IT^ be the set of increasing sequences io = < i\ < ■ ■ ■ < i p = d. 



For 



1 < i < define the partition vectors of s G C d from the sequence (ii, • • • ,i p ) to be the vectors 
s (j) := (s l3 _ 1+1 ,--- ,5^), 1 < j <p. 

The following explicit formula for the renormalization of regularized MqZV is the g-analog of 
[El Theorem 3.8] which follows from [HI Theorem II. 5.1], which, in turn, is built upon the idea of 
[II Theorem 4]. 

Proposition 4.3. Let P be the operator sending a Laurent series to its pole part: P(53 n >— JV a « e ™) = 
^2Z=-N a « e ™ an d ^ P = — P- Then for any s £ Z d and re (R>o) rf 



«:-]) 



p 



(«!, 



■■ ,i P )Gn d 



P 



(4( 



s -<2) 
f(2) 



P z 



rs <l) 

f(l) 



(33) 



We are going to use the map Z q+ to define the renormalization of MqZVs. But before doing so 
we recall that at the beginning of section [3] we mention that the multiple variable cases are different 
from the single variable case, just like the situation where a function of two variables can have all 
the directional derivatives at some point but yet is not differentiable there. Such phenomenon won't 
happen to functions of one variable. The renormalization process is essentially a limit process just 
like taking the derivatives. The behavior of multiple zeta functions at poles arc not so bad in that 
if we take appropriate "paths" to renormalize then we can produce values compatible with both 
the stuffle relations and the original values if they are defined originally. All the above remarks are 
still valid for the g-analogues. The appropriate "paths" in our case is given by Definition 14.71 and 
Definition 14. 8 1 lat er in this section. We first define the directional version of these as follows: 



Definition 4.4. For s £ Z d and r e (I 



■>0) 

s 
r 



the renormalized directional MqZV is defined by 



= lim Z„ 



and r is called the directional vector. 

Corollary 4.5. The renormalized directional MqZVs satisfy the q-stuffle relations 



^aawi?]) =<«([?]•«[?]) -£<'-«>• e<.( 

m>0 



(34) 



where 



E n >o,m>o(^ 



Proof. It follows directly from Definition 14.41 by Proposition 14.1 



□ 



The following proposition is the g-analogue of jT2l Corollary 4.13]. This and the length three 
case provide us a hint at the general shapes of the renormalization of MqZVs at non-positive integers 
which will be given in Proposition 15.41 



Proposition 4.6. For si, S2 < 0, r\, r^ > 0, set t = 1 — s% — S2- Then 



c,( 



Sl, s 2 

n,r 2 



In q 



-si 

E 

3=0 



tqjsi +S 2 -1) 
Sl -1 



-si 



1-j-Sl 

E 

i=0 



E 

3=0 

1-3- 



Sl 

3 



l=t 



l + l 

I - t 



i-t 



n + r 2 



l + si+j 



si 



(l- q )%( S2 -j-iX g (j + Sl -l). 
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Proof. From equation (|3Tj) in Corollary |3.4l we have 



Z 



■ (\ SllS2 ]\ 



E 

3=0 



-Si 



1 + 1 



ln <7 fb^ i ! 



E 

3=0 

oo ;+i 

E E 

l=—s\—j i—0 



z„ 



(!-«)%( 



s 2 -j - i 
r\ + r 2 



C g (-/)(ne) 



l+si+j 



(l + 8l+j)\ 



By Theorem [321 we get 



Si, S 2 

ri,r 2 



In 5 



= fg-iVvHi)7-iy- 



E 

3=0 



j! \ri£/ 



(n + r 2 )e 



3-S2+1 



(i-« 2 )! 



9 _1^(_£ 1 )1/_1n1-.i-j ~ . (( ri + r 2 )e) 



+ 



g-i 

lng 



E 

3=0 



■Si 



"Si 



00 Z+l 

E E 

/ — — si —3 i—0 



fc=0 

i 



iky. 



(35) 



{-(n+r 2 )ey+ 



i—s 2 + l 



(l + Sl+j)\ 



00 l + l 

E E 

i=— si —j i—0 



1 + 1 



+ E 

3=0 

By Definition 14.41 and Proposition 14. 31 we have 



fc=0 



((ri+r 2 )e) k 
(k)l 



(1-9) 



C,(-0(rie) 



i+si+3 



c 9 ([ Sl ' S2 l) 

Uri,r 2 J/ 



lim 



P z 



Si, S2 

ri,r 2 



It follows from Theorem 13.21 and the above expression for Z q ( 



Si,s 2 

T\,Ti 



e ] that 



c( 



si, s 2 " 

9-1 
In g 



g-1 
lng 



»'2\ 1 - sl Cg(si +S 2 - 1) 



E 

3=0 



-si 



E 

1=1— SI —S2 



1-si 

/ + 1 
i + Si + s 2 — 1 



3=0 

(1 - g ) , +«+«»- 



(-si)! fr-L+r 2 \ x -*i -J Cg (si + s 2 - 1 ) 
(1-si-j)! 



E l-siJ! m + r 2 \ 



\ n + r 2 / 



E 

3=0 



si / \ 1-3-Sl 
Si 



E 

i=0 



1 - j - Si 



n + r 2 - 

(l-grC g (s 2 -j-i)Cg(j + si-l). 



J + si + j 



A simple combinatorial formula quickly reduces this to the formula in the proposition. 

We are now ready to define the renormalization of MgZVs. 
Definition 4.7. For s£ (Z >0 ) rf U (Z< ) rf define 



□ 



Cq(s) = lim Cg( 



s 

\s\ + s 



(36) 



where, for s — (s±, ■■ ■ , Sd) and 5 G K>o, we write |s| = (|si|, • ■ • , \sd\) and |s| + S = (|si| + 
(5, ■ • • , \sd \ + S). These values ( q (s) are called the renormalized MqZV of the (i-tuple g-zeta function 
Cg(«l, • • ■ , Ud) at s. 

In order to deal with the g-stufne relations we need the shifted version of the about definition. 

Definition 4.8. Let / = (/1, ...,//) be a binary vector with entries equal to cither or 1. Let 
s G (Z<o)' or s G (Z>o)' such that \s\ — f G Z' >0 (i.e., if Sj = then = 0). We call such a binary 
vector / a shifting vector of s. We define the shifted renormalization of Cg(s) 6y / as the limit 



Hm ft, ( 



(37) 
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For example, ( { q hm (~l, -2, 0) = hm C, 



1,-2,0 
25, 2 + <5, 5 



where 2(5 appears because of the stuffings 



at the first position which is hinted by the shifting vector (1,0,0). 



We will compute the limit in Theorem l4.12l and show that the limits ([57)1 exist for s G (Z>o) 
in Theorem 14.131 for s € (Z<o) d in Theorem 15.21 First we consider the case of Riemann g-zeta 
function. 



Theorem 4.9. Let 7 be Euler's gamma constant and define M(q) by (|23[) . Then 

C,(i) 



In (7 



In g 



and, /or integers s > 1, Cg( s ) * s ^ e uswaZ Riemann q-zeta value ( q {s) defined by the series |3]). if 
s = —I is a non-positive integer then 



1 



t q (-l) = ( q (-l) = (l-q)- l {J2(- 1 ) r 



1 



,Z+l-r 



1 (Z + l)lng 



Proof. The expression of Cj( — is given by [TBI (60)]. The rest follows from Theorem 13. 21 



(38) 



□ 



Our primary goals are to show that our definition of renormalizations of MgZVs are well-define, 
these values agree with the usual MgZVs whenever the usual values arc defined, they satisfy the 
g-stufflc relations, and they become renormalizations of MZVs when q f 1. 

Theorem 4.10. Let s = (si,-- - , s<j) G Z> such that si > 1. Let f G (Z > o) d oe an arbitrary 
vector. Then _ J = Cg(^) which is independent of the choice of f. Here C, q {s) denotes the usual 
definition of multiple q-zeta values. In particular, Cg(s) = Qq{s) satisfy the q-stuffle relation. 



Proof. By definition Z q 



e J converges uniformly for e G (—00, 0] and therefore continuous as a 



function of e. In particular it is a power series and lim £ ^ Z q ^ * ; e 

Now we consider the divergent case s\ = 1. For f(e),g(e) G C[T][e], denote /(e) = 5(e) + O(e) 
if.g(e)-/(e)GeC[T][e]. 

Lemma 4.11. For c > set X = ^-(lnc + T). Let s G (Z> ) ! where either I = or Si > 1. Lef 
r G (K>o) z and l d = (1, 1, • • • , 1) G Z d . T/ien 



Id, s 
ci*d,r 



= P d j{X) + O(e), 



(39) 



where P c i g{X) is a degree d polynomial in X with leading coefficient C, q (s)/d\ (by convention, if 
I = we set Cq(s) = lj- Moreover, Pd^(X) is independent off. 

Proof. We prove the lemma by induction similar to that of [121 Lemma 4.4]. Notice that the terms 
produced by shifting will produce polynomials of degrees less than that of the leading term. 

When d = the lemma follows from the proof of Theorem [4.101 as eC[e] C eC[T][e]. We now 
fix d = 1 and prove the lemma by induction on I. When / = the lemma readily follows from (1271) . 
Assume equation (f3l?f is true when the length of the vector s is < / for / > 0. Then by (|2"T|) and 
Theorem 14. 101 we have 



(X + 0(e))Z q 



-Z n 



l,s 

c, f 



= u\ Z, 



1-1 



3=1 



,rj,c,r J+1 ,- ■ ■ ,n 



EM 



■z q 

ce-j 



-(1-9)EM 



s 

r + ce-j 



O(e) 
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where e-,- is the j-th unit vector of length I. Hence by induction assumption, 



l,s 

c, r 



= t q (s)X + 0(s). 



Assume equation (|39|) is true for every I > when the length of the vector (1, . . . , 1) is < d for 
d > 1. Then by dHJ) 



Id, 8 

cld,r 



(P d (X) + 0(e))(X + 0(e)) = f d +i(X) + O(e) 



for some polynomial fd+i(X) of degree d+1 with leading coefficient £ q (s)/dl, independent of r. On 
the other hand, by the g-stuffle relation 



i-i 



(d + 1)Z, 



9 U c i d+1 , 



Id, Si, ' ' ' , Sj,l, Sj+i, ■ ■■ ,3l 

cl d ,r 1 ,--- ,rj,c,r j+1 , ■■ ■ ,r t 



l d , s + ej 
eld, r + cej 



-(1-9) EM 



Id. « 
cl d ,r + cej 



By induction assumption 



ld+i, s 
cld+i,r 



= P d +i,?(X) + O(e) 



where Pd+i,s(-X") is some polynomial of degree d+1 with leading coefficient ( q (s)/(d + 1)!, inde- 
pendent of r. This completes the proof of the lemma. □ 

Now we can show 
Theorem 4.12. For sE (Z >0 )"\ ( q (s) = C g ( 1 ) ■ 

Proof. In view of Theorem 14 . 1 01 we only need to consider the case Si — 1. Assume then s = (1^, s*) 
where s* = (si, . . . , s;) and si > 1. Observe that for any substring s" of s and substring f" of r by 

the above lemma Z q [ _ 4I ) has no pole part. This means that in equation (|33J) of Proposition 
there is only one non-trivial term which gives the finite part 

Z? +([r*]) = (/rf ~ P) (M[rM) = P d,3>(X)+0(e). 



Hence 



By Definition 14771 



<.([3) 



Pd.s'(X). 

(9-1)2"- 



In q 



□ 



However, when si = 1 the g-stuffle relation are not exactly preserved under renormalization 
because of shifting. For example, Proposition 14.11 implies that Z q+ is an algebra homomorphism 
■Kz — >C[T]{e\ with 



Z, 



9+ 



Thus for c > we have 



T- 


)'-«,( 


T- 


*? 


T- 


) 


-C- 




. C- 




. c_ 





) ^z q ( 



=2Z n 



c, c 

1,1 

c, c 



e=0 



Z n 



z„ 



= C 



r 2 • 

2c. 
2 

2c. 



+ (1 - ?)^ 
+ (1 - q)Z q 



(l-g) 2 ln2 
lng 
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by Theorem [372] This implies Cg 1} (l) = CgCQ + (1 - g)ln2/lng and 

=2C g (l, 1) + C g (2) + (1 - g)C g (l) 
C 9 (l) 2 =2C g (l,l)+C ? (2) + (l- 9 XW(l). 
Of course, when q f 1 these two values are both equal to C(l) 2 = T 2 . 

Theorem 4.13. Let s G (Z >0 ) i with S\ > 1, then for any positive integer d the limit of CqO-d,s) 
in (|37l) exists. Moreover all the shifted renormalization of MqZVs at all positive arguments satisfy 
the q-stuffte relations in the following sense. Let s% £ (Z>o) fc , s*2 £ (Z>o)^ and assume there are Wt 
vectors of length k + I — t produced by the q-stuffle with t stuffings and let fj be the binary vector 
representing the stuffing positions of the j-th such vector: 



min(fc,f) W t 

Sl* q S 2 = E (l-Q^E^ 
t=0 3=1 



Then 



m'mk,£ Wt 

C,(Si)Cr(Sb) = £ (i- E $ - /I 

i=i 



4=0 

Proof. By the very definition of the shifted renormalization and the q-stuffle relation it suffices to 
show the first part of the theorem, namely the existence of the limit in (|37|) . Indeed, if the leading 
component in s is greater than 1 then £ 9 is independent of the directional vector by Theorem l4.10l 
In the following we consider the shifted renormalization Cg only. Let the vector r = 

(r%, ... , rk+i) = (la, s) + / be the shifted directional vector of length I + d. Then 



Id, s 
r 



e n 

kl>»->k d j=l 



exp(r j £[fc 3 ]/g fc 



e n 

fc (J >ni>-->n ! j'=l 



!) exp(r J £[nj]/g™ ; ') 



(40) 



We want to show this series is good (used just in this proof) in the sense that it is in C[e][ln(— e)] 
and is finite when n, • • • , are equal to either 1 or 2. We prove this by induction on d. If d — 1 this 
follows from (|3"5|) in Lemma H. Ill Assume the series in (|4T)]) is good when the length of (1, . . . , 1) 
in front of s is less than d for some d > 2. Set as before 



F(x,s;e) 

As in the proof of Theorem 
expression for 



jc(s i) ex p(eq x [x\) q x( - s 1%> exp(e(q x — 1)/(1 — g)) 



(l-q) s [x] s (l-g x ) s 

we can use Euler-Maclaurin summation formula to find a closed 



EnM;e) = E 



fc>i 



exp(e[/c]/(7' c ) 



111 

F(as + i, 1; e) dx + -F(i + 1, 1; e) - —F'(i + 1, 1; e) - - 



The following integral 



+ i, 1; e) da: 



can be evaluated using (j2"l~j) by the substitution e — » e[z + l]/g l . Therefore we get 



B 2 (x)F"(x + i,l;e)dx. (41) 



(1-9) 



F(x + i, 1; e) dx = 
1-9 



1-9 



In g 



lng 

ln(— e) — In q + In 



In 



-e[i + l] 



+ 7 + E 



iu q i(i+i) 

i 00 / 

h ( 1 + M w ) + ^gw'( 



[i] 



(42) 
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Let i = &2 in the above two formulas (|41j) and (|42j) . It is straightforward to see that the middle two 
terms of (|4T|) both contribute to good series by Theorem 14. 101 and the induction assumption. The 
last term of (|4Tj) can be handled similarly using (|16|) after we notice that P>2(x) is bounded by 1/6 
from an easy computation from the series expansion (|20|) and the fact that £(2) = 7r 2 /6. So the 
main divergence term of (|4"tl)) when e is near zero comes from the first term of (|4"Tj) . Thus by (|4"!2)) 
it is enough to show that all the sums below are good: 

z ,„(M)n "-<^y/«") £ (43) 



k 2 >--->k d ' j=2 JJ <!d>iii>->iii 



gmfc2 II eXP(rj ff /gfc3) E -,™^o (44) 



^ \k 2 ] m 11 fc, 

k 2 >--->k d L J 3=2 L JJ fc<j>«i>— >n 

e'[fe] m i — r exp(rje[kj}/q k 3) 



H ^jwr; ... )1<m</ . ( 45) 

k 2 >--->k d ^ 3=2 L JJ fe<i>ni>— >ni 

In if m = then this follows from induction and if m > then it follows from Theorem 14. 101 
For (|45p . we can mimic the argument for finding Z q (s, e) at non-positive integers by differentiating 
and see immediately that ([43]) is actually in C[e] and r 2 can appear in the denominator only 
in the form of some pure power. So the series in (|45[) is good by induction. 

The most difficult one is (j43[) in which case we again apply Euler-Maclaurin summation formula 
using the following modified version of F(x, s; e): 



qx{s l) ex p( £ g x [x\) 
{l-q) s [x] s ' 



Then ([42]) changes to 



G?j-(x + i,l;e)da;= / ( ^ (lni) J i~ 1 e~ t rft 



1 , ln fe+i f £ [* + 1 h__/ ,e [* + 1 h , 1 



j + 1 V / V <f / j 



(46) 



This integral can be treated similarly as (|42|) by using polygamma functions (see [T]) which is closely 
related to 

/■OO 

f k (s) := / (lnt) k t 1 - s e- t dt = T^(s). 



JO 

In particular, to proceed by induction we need the fact that all the values of /fe(l) are finite which 
are in fact Q-linear combinations of the weight k products of C(n)'s an d the Euler constant 7, 
where we take the weight of 7 to be 1. For example, the first few values are /i(l) = —7, /2(1) = 
C(2) + 7 2 , / 3 (1) = -2C(3) - 3 7 C(2) - 7 3 , and / 4 (1) = f C(4) + 8 7 C(3) + 6 7 2 C(2) + 7 4 - This shows 
that the sum in (|43[) lies in C[e] and the only form that r^'s can occur in the denominator of any 
coefficient of the series is in a factor of some form r it + • • • + r it . Note that rj's can also occur in 
logarithms in the form ln(rj) which implies that we can take r^'s to be either 1 or 2. This finishes 
the proof of the theorem. □ 

5 Renormalization of MgZVs at nonpositive integers 

The following result is straightforward: 
Theorem 5.1. For s € (Z<o) d , we have 

M = C([J r J)= ( £m^[3). 
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Proof. This follows from Corollary 13. 5i Definition 14.41 and Definition 14.71 

In the rest of the section we are going to prove 
Theorem 5.2. Let s £ (Z<o) d and f be a shifting vector of s. Then the limit in (|37[) 



□ 



lira Cq( 



S+\-s+(S-l)f\ 



exists. 



The major complication of the proof of Theorem 1 5 . 21 arises from the possibility that Sj can be 
or — 1. To prove the theorem we need some more information of Z q ^ 



s 
f 

Definition 5.3. Let r £ (M > o) <i . Let d and j be two positive integers such that j < d. Let 

N(f) =rt+ r i+ i + - + r d , D{f) = r,- + r j+1 H \- rj, for 1 < j < i < d. Then we say N(r) w /D(r) v 

is a tailored fraction for any nonnegative integers v and w such that w > v. For convenience, we 
regard all polynomials of r as tailored fractions too. If a rational function P/Q is a product of 
tailored fractions as above then we call it a tailored rational function. Note that deg(P/Q) > 
always holds. If every coefficient of a Lauren series Yl c i £% IS an R-lmear combination of tailored 
rational functions then we say the Lauren series is a tailored Laurent series. 

Proposition 5.4. Suppose m and n are two integers and d is a positive integer. Let s £ (1<o) d 
and r £ (K>o) d . For (i\, ■ ■ ■ ,i p ) £ lid, let 1 < j < p be the partition vectors of (s2, ■ ■ ■ , sd) 
from (ii, ■ ■ ■ , i p ) in Deflnition \4:.2\ and similarly define ^\ 1 < j < p. Let 



i>o (ii,- ,i P )en d 



r gip) 
f<p) 



■P z, 



rs <2) 
^2) 



P e m Z 



rs <D 



(47) 

Then for every i the coefficient Ci can be expressed as an M.-linear combination of rational functions 
P/Q o/ri, • • ■ , rd with P and Q having no common factors such that the following must hold: 

(a) P and Q are homogeneous polynomials such that deg(P/Q) + n + m = i. 

(b) Either Q is a constant or every factor of Q has the form rj + ■ ■ ■ + Td for some j = 1, . . . , d. 
So Q is uniquely determined if we require Q to be a monic polynomial with respect to r d . Then 
either Q = 1 or Q is a product of linear factors of the form rj + ■ • • + r d (j — 1, . . . , d) with possible 
repetitions. 

(c) If m = then rd divides Q only if n > 1 — Sd- 

(d) If m = and n < then P/Q is tailored. 



rg <t) 



5^ ai£ l is a homo- 



Proof. (a) By Corollary 13.51 for 1 < t < p every coefficient a, of Z q (^ 

geneous rational function P/Q of r\, ■ ■ ■ ,ra such that deg(P/Q) + i = 0. Now the product of two 
such terms 



{Px/Qxy ■ (P 2 /Q 2 )e J = (P 3 /Q 3 )e 

still satisfies deg(P3/Q3 + i + j = 0. Part (a) follows easily. 

(b) We proceed by induction on i If d = 1 then this follows immediately from Theorem 
Now assume part (b) is true when the length of the vector s is d— 1 for some d > 2. Let s £ (Z<o 
We have the disjoint union 



n d ={(i,ii + i,. • • 
u{(ii + l,i 2 - 



i p + 1) | (h, ■ 



(h, 



g n d _i} 

• •• ,i p ) £ n d _i}. 



Note that i p = d—1 in the above by definition. In the rest of the proof of (b), for («!,••• ,i p ) £ ILj_i 



let 1 < j < p be the partition vectors of (s2, . . . , Sd) from (i\ + 1, • • 



1) in Definition 1421 
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and similarly define , 1 < j ' < p. Putting = ■■■ i )en d x we nave 



rs -<P) 



r|(l) 



r 

By Theorem l3.2l we have 



p 



(>%( 



o 

g -iA(n)(-«i)l 



where the sum is vacuous if m > and 

A(r l): =-(^) 
Write fc = ii + 1. Then by (J3TJ) we get 
g-lB(rx)(- Sl )! 



E CrC-i-^ 

j=0 



* 1 1 



if TO > —si, 
if TO < —Si, 



1\ 1-si 



(48) 



lng e i-(m+ Sl ) 



E 



-si 



E § i + % 



S2 -J,S 3 ,...,S k 

r\ + r 2 ,r 3 , . . .,r k 



(l + si +j)le~ m 



where 



^):=E^i 



3=0 



Si -j,s 3 ,...,s k 
r\ +r 2 ,r 3 ,.. . ,r k 



e\e J . 



(49) 



By induction assumption, the linear factors of Q can only be r\, or r,j + • • • + r c [ (j = 1,2, ... ,d). 
Let us exclude the possibility of r\ . It suffices to prove that 



lim A{ n )Z q 

ri— »0 1 



^1) 



Set 



m=z q ( 



S 2 , S 3 , ... , Sfc 

x + r 2 ,r 3 ,...,r k 



- B(r\) J is finite, 
ej . Put D = d/dx and £ = — s% > 0. By L'Hopital's rule 

/(o)-£i(-*)^7(z) 



lim /(O)A(x) + B(x) = (-1)* lim ■ 

x^O x— >0 



j=0 ' 



pt+1 



-l) 4 lim— 



(D t+1 f)(0) 

(t+iy. 



< oo 



after cancellations in the numerator. 



Remark 5.5. By exactly the same argument we see that r<2 can not appear in the denominator of 
flggP when d > 2. 

(c) and (d). Assume to = 0. We use induction on d again. Theorem 13.21 yields the case d = 1 
easily for both (c) and (d). Suppose part (d) holds when s has length d — 1 for some d > 2. Now 
let s £ (Z<o) d . Observe that we have another disjoint union of given by 

n<j ={(ii,--- ,i p ,d) | (£].,••■ ,i p ) € n d _i} 

U{(*i,--- + 1) | • • • ,*p) € n d _!}. 
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Similar to part (b), in the rest of the proof, for any positive integer t and • • • ,i p ) € Tit, let 
1 < j < P be the partition vectors of (si, . . . , St) from • • • , i p ) in Definition 14.21 Define r^' 
in the same fashion. Put £ r(() = ,» p )en t ' 



Then 



r(d-i) 



r3(l) 



E 



P 



s -<P-i) 



•P U 



For ease of notation, in the rest of the proof we write P,; = r,; + • • • + r c i, k = i p -\ + 1, and 



X = P 



(«,( 



rg G>-i) 



•P U, 



/7-n.o/ 
Z 9+ \ 



i(d-l) 



Using (f3"2"]l we get 

S<P) 



r(d-i) 



lno ^7! ^ \ Vr^e/ 

i 3=0 J T(d-l) 



Z„ 



Sfc, • • • , Sci-2, Sd-1 - J 

^fcj •• • ,r<i-2, Rd-i 



e\X 



E 

3=0 



-Sd 

j 



E E^ 

(=— Sd— j i(d—l) 



,n<q{~l){rd£) 



l+Sd+3 



(1 + Sd+jY- 



ol+l 7 



Sfc, ■ • ■ , Sd-2, S d -1 - J 
rk: ■ ■ ■ ,Td-2, Rd-1 



(50) 
(51) 
(52) 
e)X). (53) 



By Theorem 13.21 and Remark 1 5. 5 1 we see that can appear essentially in the denominator of some 
term in (|50|) only if n > 1 — Sd which proves (c). As a matter of fact, we don't need induction 
assumption to prove (c) so we will use (c) freely in what follows. 

Suppose now n < 0. By part (b) every denominator in (|51j) can only have linear factors of the 
form Td or Rj — r& for some j = 1, . . . , d — 1. By part (b) and (c) we know that these terms will 
be cancelled out in the end so (|5ip doesn't contribute to any un-tailored terms. In fact, negative 
powers of Rj — Td can not really appear by the reasoning in the next paragraph. 

As the exponent of e in front of Z q is now negative in (|52| every term has Td in the denominator 
so that all of these will be cancelled by corresponding terms in ([52]) (note that there is no negative 
powers of in (|53[) ) . In fact by part (b) the only linear factors of denominators that can appear 
are of the form Rj in ([52]) and of the form either Rj — rd or in ([ST]) . So to cancel terms with 
negative r,j-power no factors Rj — r& can appear in denominators of (|51[) and no factors Rj can 
appear in denominators of (|52j) . (Actually, positive powers Rj and Rj — Td must appear so that 
Td = Rj — (Rj — Td) is produced on the numerator to cancel negative powers of Td)- This shows 
that the sum of (fBTj) and (|5^|) is tailored. 

Let's turn to (|53p . If I + j is small, say I + j < —Sd — n then by induction assumption the 
Laurent series is tailored. When I + j increase gradually it seems that un-tailored terms may 
appear. We will show this actually can not happen. Roughly speaking, in order for some linear 
factor rj + r^+i + • • • + r<j to appear in the denominator, the power of e in front of Z q in (|53|) has 
to be large by part (c). But this will produce high powers of r<j which results in tailored fractions. 

It is easy to see that in (|53l) the highest power of Td-i + Td that can appear in the denominator 
is 1 — Sd-i + j by part (c). This takes place if I + Sd + j + n > 1 — Sd-i + j- Similarly, we can 
expand (|53|) by iteratively using the formulas from (|50|) to (f53|) and see that the highest power of 
Ri that can appear in the denominator is 1 — Sj + jd-i where jd-i is the corresponding summation 
index in ([52]) and ([53]) . Indeed iterating once we see that for fixed j and I the essential part of ([53]) 
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becomes: 



l+Sd+j 
' d 



P(e n+l+Sd+3 Z q ( 



i(d-l) 



Sk, 

Ofe, 



, S d -2, Sd-l - J 
■ , Td-2, Rd-1 



r d +Sd+3 P(e n+l+Sd+ JZ q 

i(d-2) 



(\ s d-i — i] \ 
U R d -i \) 



\P z, 



s)X 



X 



1 3Sd-l , . 



lng 



J2=0 



J2! 



i*(d-2) 



(54) 
(55) 

(56) 



E 

J2=0 



oo 

E E 



1 \ l+i-sd-i— n 



Sk, ■ 

rk, 



, rd_3, i?d-2 



J - Sd-l 

h 



(^2 + S<J - i+j2)! 



h=j— Sd-l— J2 2) 



(57) 



Sfc, ■ 

rk, 



, Sd-3, — 32 

,rd-3, Rd-2 



e)X 



For (|SSJ) observe that 



|(P) 
ftp) 



)Xj can not produce legitimate un-tailored fractions by part 
(b). So the only possible un-tailored terms are the result of negative powers of Rd-i coming 



from Z„ 



Rd-l 



When n < the only way Rd-i can appears in the denominator is when 



n + l + Sd+j>l — Sd-i- Then these terms are multiplied by rJ~ Sd+J to become tailored. 

Let's consider (|5T))) next. Put w := n + 1 + Sd + j — (1 + j — Sd-% — 32)- We treat the two cases 
w < and w > separately. If w < 0, then the induction assumption takes care of all the terms 
except negative powers Rd-i- But by breaking ([54"]) into two parts like in ([50)) we see that when 
w < negative powers of i£d— 1 can n °t appear by part (c). So we can assume now w > 0. As n < 
this means that l + Sd + j>n + l + Sd + j>l+j — Sd—i — 32 which implies that 



J+Sd+j n+l+s„ 

1 w 



\Rd-iz) 



1 \1+J-Sd-1-J2 



Td 
Rd- 



l + Sd+j 



£ n (Rd-ie) 1 



is tailored. Hence we see that (|56| is now tailored by setting a = w — n > and £ = d — 2 in the 
following Claim. 

Claim. Let a and t be two positive integers such that 1 < t < d — 1. Then the series 



^P(>(i* t+ie )<%( 



Sk, 

rk, 



■ , St—l, St 

■ ,r t -i,R t 



s)X 



is tailored. 



We use induction on t. Notice that this is the inner induction loop. The outer induction loop is 
on d. When t = 1 the claim follows from Theorem 13.21 easily. This already proves the proposition 
if d = 2. We now assume d > 2, t > 1, and the claim is true if the length of the vector in the claim 
is less than t for some 1 < t < d. 

As before we may break lit into two parts and form the disjoint union 

lit ={{h, ■ ■ ■ , i P ,t) j (h, ■ ■ ■ , i P ) G n t _i} 

+ 1) I (n, •■• ,i P ) g n t _i}. 
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Adopting the same argument we used to obtain ([50]) to (f53|) we have 



^P(>(iW)%( 



Sk, 



j St-li St 

,r t -x,Rt 



e)X 



i(t-i) 



2 p(e"(ii t+ie rz,([^])p(z g ([! (p) ])A:) + e »(J2 H . ie )«Z,( 



2 p(e"(i?, +l£ ) Q Z f; ([^])p(z 9 ( 



g-1 (-*)! 



In o ■« — ' 7 ! 

H 3=0 J T(t-1) 



E n^^H^ 



1 \ 1-at-J 



Sfe, 



,Sf_ 2 ,St-l — J 

.,r t _2, -Rt-i 



E 



i 



oo 

£ 2 P(e n (i^iE) 



tt C 9 (-0(^e)' + * t+, 'g I+ i 



i=-»t-j ?(*-!) 



G + Si+i)! 



Sfe,...,S t _2,S t _i — J 

r fe ,...,r t _ 2 , i? t -i 



(58) 
(59) 
(60) 
(61) 

(62) 



By part (c) i?t appears in some denominator of the first term of f|59[) only if n + a > 1 — s t in 
which case the exponent of Rt is l — s t . Then multiplied by this becomes tailored as a > a + n. 
Un-tailored terms with R t appearing in the denominator can not be produced by the second term 
by Remark 15.51 

The same argument for (l5Tj) above using part (b) and part (c) shows that (|60p does not contribute 
to any un-tailored term in the end except possible powers of 1/Rt- However, all such terms are 
multiplied by high enough powers of Rt+i resulting in tailored terms again. 

Similar argument for (|52p implies that in (|6ip we only need to consider the case when n + a > 
1 — St- In this case we see that 

-1 \i-*t-j 



i?i+l\ a 
i 

is tailored and therefore (|61|) is tailored Laurent series by inner loop induction assumption. The 
same argument by manipulating the powers of e works for (|62|) in the same fashion and this finishes 
the proof of the claim. 

Finally, ([57)1 follows from the claim immediately by setting a = l + ^s^ + l^+ji > and t = d— 1 
since ^ 

\Rd-lJ 

This finishes the proof of the proposition. □ 



Proof of Theorem 15.21 Clearly Theorem 15.21 is true when d = 1 by Theorem 13.21 When d > 2 by 
setting m = n = in Proposition 15.41 in (d) we see that the constant term cq of Z®+^ ^ must 
be the product of linear tailored fractions of the form 

n + r i+1 j h r d 

r i + r i+i H 1" 7 'd 

for some i and j such that d > i > j > 1. Note that in the limit ([3~7[) each can only be either 
5 + rik, or 2(5 + rtfe for some nonnegative integer rik ■ If nk > for some j < k < d then the limit 

,. r. + r, + H hr d 

lim 



<5^o+ Tj + r i+ i H h r d 



is clearly finite. If rik = for all fe = j, j + 1, . . . , d then the limit is still finite because S will be 
cancelled out. 

This completes the proof of Theorem 15.21 □ 
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One of the most important properties of MZVs and MqZVs is that they satisfy the (g-)stume 
relations. In order to show that our renormalization for MgZVs at all non-positive arguments is 
correct we need to show some type of q-stuffle relation hold for them. However, we must use the 
shifted renormalization following the shifting principle of MqZVs. 

Theorem 5.6. The shifted renormalizations of Cq{s) satisfy the q-stuffte relations for all s£ Z< . 
Explicitly, let s*i £ Z< , 5*2 £ ^<o an d assume there are Wt vectors of length k + £ — t produced by 
the q-stuffle and let fj be the binary vector representing the positions of the stuffing in j-th such 
vector: 

min(fe,£) W t 

Sl* q S 2 = ( X ~ J2(Vj - fj). 

t=0 j=l 

Then 

minfc,^ Wt 
t=0 j=l 

When q] \ we obtain \1'A Theorem 4-11]. 

Proof. By the definition of renormalization in Definition 14.71 and Corollary 14.51 we have 
C q (sx)C q (s 2 )= lim ^([l 1 ])^! 2 !) = Jim c( 







"s 2 - 








-f\- 




■r 2 - 



From the definition of q-stufflc (flO| and the shifted renormalization (|3T[) the theorem follows im- 
mediately. □ 

To conclude our paper we observe that a multiple q-zeta function has more singularities than 
its classical counterpart. We don't know how to renormalize MqZVs at these points at present. 
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